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Abstract 

Low energy effects induced by heavy extra degrees of freedom are suppressed by 
powers of the large mass scale, thus preserving, if sufficiently heavy, the successes 
of the Standard Model in describing low energy phenomena. However, as is well 
known, heavy right-handed neutrinos may play an important role in low energy 
,^ I phenomenology as an explanation of the smallness of neutrino masses. We consider 

in this paper an extension of the Standard Model by heavy right-handed neutrinos 
and heavy Higgs doublets and we show, using a renormalization group approach, 
^ i that this model can simultaneously provide an explanation for the small neutrino 

QQ I masses and for the mild hierarchy observed between the atmospheric and the solar 

cn ■ mass splittings, even when the extra degrees of freedom are very heavy. We analyze 

the necessary conditions to reproduce the oscillation parameters and we discuss 
t^ I possible experimental signatures of this model. 



X 



1 Introduction 



The last fifteen years have witnessed a tremendous experimental progress in neutrino 
5^ \ physics, leading to a good determination of two mass splittings, two mixing angles and 

a fairly stringent bound on the third angle [1]. The most conspicuous difference between 
quark and neutrino properties is without any doubt the smallness of neutrino masses. 
Furthermore, the better and better measurements of neutrino parameters have also re- 
vealed the existence of large mixing angles in the leptonic sector and the existence of a 
relatively mild mass hierarchy between the two heaviest neutrino masses [2]. 

Extending the particle content of the Standard Model (SM) with three heavy right- 
handed neutrinos, thus implementing the (type I) see-saw mechanism [3], solves very 
elegantly the problem of generating small neutrino masses and opens new opportunities to 
understand the puzzles of the existence of large mixing angles and a small mass hierarchy. 
Unfortunately, whereas the seesaw mechanism does not make any generic prediction 



about the leptonic mixing matrix, it tends to predict a neutrino mass hierarchy which 
is much larger than the one inferred from experiments [1]. Namely, under the plausible 
assumption that the neutrino Yukawa couplings have hierarchical eigenvalues, as observed 
in the quark and the charged lepton sectors, the mild mass hierarchy observed between the 
solar and the atmospheric mass splitting can only be accommodated in very special cases. 
One possibility arises when the hierarchies between the masses of the heavy right-handed 
neutrinos is much larger than the hierarchy between the neutrino Yukawa eigenvalues, 
in which case the right-handed mixing angles have to be tiny. Alternatively, when the 
hierarchy in the masses of the heavy right-handed neutrinos is comparable to the hierarchy 
in the Yukawa couplings, it is possible to accommodate the observed mild neutrino mass 
hierarchy only for certain, possibly fine-tuned, choices of the right-handed mixing angles. 

Further extensions of the type I see-saw model have been considered in the literature. 
A minimal possibility consists on introducing one extra Higgs doublet, which leads to 
radiatively induced neutrino masses. Most works impose an additional Z2 symmetry, 
under which the Standard Model particles are even, whereas the right-handed neutrinos 
and the extra Higgs doublet are odd [5l|6]. With this assignment, the tree level neutrino 
mass vanishes and the leading contribution is hence the radiatively generated one. In this 
scenario, the right-handed neutrinos and the extra Higgses could be directly produced 
at colliders while being the neutrino masses in the measured range, thanks to the loop 
suppression and an appropriate choice of the neutrino Yukawa and Higgs self-couplings. 
Furthermore, the Z2 symmetry ensures the absence of tree level flavour changing neutral 
currents. On the other hand, it has also been discussed the general two Higgs doublet 
model (2HDM), without imposing ad-hoc discrete symmetries |^. In this case, both 
mechanisms of neutrino mass generation are present, naturally leading to a mild neutrino 
mass hierarchy [S]. 

In this paper we carefully analyze the mechanism of neutrino mass generation in 
the two Higgs doublet model extended with right-handed neutrinos. Rather than being 
motivated by finding signatures of new physics in experiments at the energy frontier or at 
the intensity frontier, we are motivated by constructing a simple and natural framework 
capable to explain the observed neutrino parameters while preserving the successes of the 
Standard Model. This approach is inspired in the high scale see-saw mechanism, which 
despite its well known lack of testability, still stands as the most compelling explanation 
for the small neutrino masses. 

We will argue that by making all the new particles heavy it is possible to simulta- 
neously explain the smallness of neutrino masses and the mildness of the neutrino mass 
hierarchy, without jeopardizing any of the successes of the Standard Model. Remark- 
ably, in this model only one right-handed neutrino suffices to generate two neutrino mass 
scales: the atmospheric neutrino mass scale will be generated at tree level, whereas the 
solar mass scale will be generated by the radiative corrections to the effective neutrino 
mass matrix. Both neutrino masses are suppressed by the large right-handed neutrino 
mass scale, thus explaining the tininess of neutrino masses. On the other hand, the ratio 
of the two mass scales is suppressed by the loop factor and enhanced by a large logarithm 



of the ratio of the heavy right-handed neutrino mass to the Higgs mass, resulting in a 
neutrino mass ratio which can be roughly of the correct size. Finally, the decoupling 
of the extra Higgs degrees of freedom ensures the absence of large contributions to the 
flavour and CP violating processes, both in the leptonic sector and in the quark sector, 
as well as to the electroweak precision measurements. 

In Section [2] we review the basic features of the two Higgs doublet model, the various 
problems which arise in this very minimal extension of the Standard Model, and how 
they can be circumvented altogether by decoupling the extra scalar degrees of freedom. 
In Section [3] we show, using a renormalization group approach, that even in the decou- 
pling limit the extra Higgs particles can play an important role in low energy neutrino 
physics, as an explanation for the mild hierarchy observed between the atmospheric and 
the solar neutrino mass scales. The viability of this model requires at least one right- 
handed neutrino and two Higgs doublets; in Section H] we comment on the differences 
of this framework with another minimal framework of neutrino masses, namely the two 
right-handed neutrino model with just one Higgs doublet. In Section [5] we calculate the 
corrections to the leptonic mixing matrix induced by quantum effects, and we argue that 
a non-zero ^13 is generically expected, as well as a a deviation from the maximal atmo- 
spheric angle which is correlated to the angle ^13. Lastly, in Section [7] we present our 
conclusions. 

2 Benefits of the decoupling limit of the 2HDM 

We consider an extension of the SM consisting on adding to the particle content one 
additional Higgs doublet, with identical quantum numbers as the SM Higgs doublet. 
The most general Higgs potential reads [9l [10] : 

V = mii$J<I>i + m22$2'^2 - [mi2$l$2 + h.c] 

iAi($l$i)2 + 1X2(^1^2? + A3($l<I>l)($^<I>2) + A4($l$2)($^$l) 
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+ 



^A5($I$2)' + A6($I$l)($I$2) + A7($5$2)($I$2) + h.C. 
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Despite being such a simple extension of the Standard Model, the introduction of a 
second Higgs doublet in general jeopardizes many of the successes of the Standard Model. 
More concretely, the Higgs potential has now a richer structure including electrically 
charged directions, which may lead to the spontaneous breaking of the electromagnetic 
U{1) symmetry if there are minima along those directions. Besides, the extra Higgs 
doublet contributes to the oblique parameters S, T and U [11] , possibly leading to values 
in conflict with electroweak precision data. Lastly, the new, in general flavour violating, 
couplings of the fermions to the extra Higgs doublet could lead to too large flavour 
changing neutral currents or lepton flavour violation. 

It is interesting that all these problems are simultaneously solved in the decoupling 
limit of the two Higgs doublet model, which is defined as the limit where one of the 



Higgses is kept light, with a mass comparable to the Z-boson mass, while the rest acquire 
masses much larger than Mz- 

To show the absence of charge breaking minima in the decoupling limit we will work 
without loss of generality in the Higgs basis where mi2 = 0. Then, in complete analogy 
to the Standard Model, we postulate the absence of unbounded from below directions 
in the potential. Furthermore, and also in analogy to the Standard Model, we require 
that one mass squared, say m^^, is negative, to allow the spontaneous breaking of the 
electroweak symmetry. Lastly, a necessary condition for the decoupling of the second 
Higgs is ml^ > |m^J/Ai > [12]. 

To study the minima of the potential, we will exploit the SU{2)l invariance to express 
the Higgs fields as 

1 f«\ ^.^±(- 



where a and (p2 are complex fields, while (fi is a real field. 

We consider now the direction in field space (p2 = CLfi, cr = bipi, along which the 
potential reads: 

Va,bi^i) = -m^J(/?ip + -m22(|ap + IfeHlv^iP + const, x |(y9i|^ . (3) 

Clearly, for any minimum of V with vacuum expectation values (v^i), {'P2), (c) there is 
a choice of the parameters a and b such that Vafi{(pi) has a minimum at {(pi), concretely 
for a = {(P2) I ifi) and h = {a) / {(pi). This is only possible, though, if the quadratic part 
of eq. (jHI) is negative, thus for any minimum of V it must hold that: 

l(^2)P + IWP<^K,|. (4) 

Utilizing this inequality it is straightforward to determine the minimum of V. The 
differentiation with respect to (fi yields 



91777^ I 
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Ai 

Varying now V with respect to a and ip2, it can be checked that 

(a) = , 



(5) 






2m^2 



which shows that the electric charge is conserved. Furthermore, as the decoupling limit 
is approached, the Standard Model vacuum is recovered. 



Apart from recovering the nice features of the Standard Model vacuum, when taking 
the decouphng hmit also the successful predictions of the Standard Model in electroweak 
observables and flavour physics are recovered. Namely, the existence of an additional 
Higgs doublet introduces contributions to the oblique parameters S, T and U which can 
be in tension with the electroweak precision measurements. However, it can be shown 
that in the decoupling limit S, T and U scale with |"^ii|/(Aim22) [13], thus bringing the 
oblique parameters within their experimentally allowed values for sufficiently heavy extra 
scalar states. 

Besides, the general 2HDM induces in general too large rates for the flavour changing 
neutral currents and the lepton flavour violating processes. Whereas this problem can be 
alleviated by assuming concrete flavour structures of the Yukawa couplings [H] , a simpler 
way to suppress altogether all new contributions to the flavour violating processes consists 
on assuming that the new scalar particles are all very heavy. 

The general flavour dependent part of the Lagrangian reads 



-C 



Yuk 



{Y^^)ijlLieRj<^a + {Y;^)ijqLiURj^a + (Y^)ijqLidRj<^a + h-c. 



(7) 



where i,j = 1,2,3 are flavour indices, a = 1,2 is a Higgs index and ^a = ^T2'^a- ^^ '^i^l 
be convenient in what follows to work in the Higgs basis where one of the Higgs fields, 
say $2; does not acquire a vacuum expectation value. In this basis, then, the Yukawa 
matrices Y^^ ^ are proportional to the fermion mass matrices. 

Consider for illustration the contribution of the second Higgs doublet to the pro- 
cess /i — >■ 67, whose rate is strongly constrained by experiments BR(/i — ;■ 67) < 1.2 x 
10"^^ [T3]. For a wide range of parameters, this process is dominated by the two loop 
Barr-Zee diagrams [TB]. The leading contribution comes from the top quark unless Y^^^ 
is small [17]. In the decoupling limit the branching ratio reads: 



BR(/i->e7) 
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Here, rrit denotes the top quark mass, f{z) is defined in [T7] and evaluates /(2) k. 1, and a 
is the Higgs mixing angle, which reads cos a ^ IW^"^ I'i^'h i^ the decoupling limit. Thus, 
the stringent experimental bound on BR(/i — )■ e 7) can be evaded, for the concrete fiavour 

if the scale of the heavy Higgs mu ^ 2 TeV. 



structure Y^^^ 



-'m33' 



^ell^e22 ^'^'^ ^«33 

Clearly, the rate of /i — )■ 67 can always be suppressed, regardless of the fiavour structure 
of the Yukawa couplings, by increasing sufficiently the mass of the extra scalar degrees 
of freedom. 

In the quark sector, strong constraints come from the measurement of the mass dif- 
ference in meson anti- meson systems, such as 5° - 5°, which arise at tree level in the 
general 2HDM. In the decoupling limit, the 5° - 5° mass difference approximately reads: 



Ami 



Am-^ + ^msJlP,--^^ 



(9) 



Where fs is the S-meson decay constant, /b^ = 238.8±9.5 MeV [18], ms, is the Bg meson 
mass, TTiBg = 5.37 GeV [19], and the coefficient Pg^^ includes the renormahzation group 
evolution from the scale Mz to ~ msg and the hadronization of the quarks to mesons [2U] 
and reads P2^ ~ 3.0. Assuming jy^sl ~ 1^^321 ~ V^d22^d33^ "^^ estimate that the 2HDM 
contribution to the B^ - B^ mass difference lies within the theoretical error of the SM 
calculation, Am|^ = (135±20) x 10"^^ GeV [21], for m^ > 3 TeV. As before, the 2HDM 
contribution to the meson-antimeson mixing can always be suppressed for a sufficiently 
large uih, regardless of the flavour structure of the quark Yukawa couplings. 

3 Neutrino masses in a 2HDM extended with right- 
handed neutrinos 

We will consider in this paper an extension of the Standard Model consisting in adding one 
extra Higgs doublet and at least one right-handed neutrino, singlet under the Standard 
Model gauge group. We will not impose any discrete symmetry on the model. Then, 
compatible with this matter content, the most general Lagrangian reads: 

£ = C^"' + C^''^ + C -V , (10) 

where £'^™ contains the kinetic terms, C^^^ is the Yukawa Lagrangian for the Standard 
Model fermions, given in eq. ([7]), l^ is the Higgs potential, given in eq. ([T]) and C is the 
part of the Lagrangian involving right-handed neutrinos, given by: 

- £" = {Y^)ijluywAa - -Muijiygi^Rj + h.c. (11) 

This term contains a Yukawa coupling, which leads to Dirac neutrino masses, and a 
Majorana mass for the right-handed neutrinos, with a size which is a priori unrelated to 
the electroweak symmetry breaking scale. 

We will assume that the mass scale of the right-handed neutrinos is much larger 
than the electroweak symmetry breaking scale and the mass of all the extra Higgs mass 
eigenstates H'^, A^, H^, which we denote collectively by rriH- Hence, the right-handed 
neutrinos are decoupled, leading to the following effective operators: 

-^'' ^'^ = +l^t(^u^a){^IlZ) + h-c. (12) 

where, at the scale of the lightest right-handed neutrino, 

k'^\M,) = {Y:M^X'^){M,). (13) 

Since we have chosen to work in the basis where ($5) = v/y/2, ($2) = 0, the neutrino 
mass matrix at the scale of the lightest right-handed neutrino depends just on the coupling 



K^': 



M.{M,) = yKi^(Mi) , (14) 



which is diagonahzed in the standard way: 

Ai„ = U*dia,g{mi,m2,m^)U'' . (15) 

These are not, however, the neutrino parameters measured by experiments, where the 
energies involved are much smaller than the right-handed Majorana mass scale. In order 
to compare the predictions of the model with low energy experiments we will make use of 
the Renormalization Group Equations (RGE) given in the Appendix, to run the effective 
couplings K,""^ from the scale Mi to the scale tuh- Below the scale rriH the neutrino mass 
matrix runs with the RGBs of the Standard Model extended with massive Majorana 
neutrinos. Since in this framework the neutrino masses are hierarchical, the running will 
not introduce any new qualitative feature but will only modify the values of the mass 
eigenvalues and the entries of the leptonic mixing matrix by a small factor, proportional 
to the tau Yukawa coupling squared and to the small logarithm \og{mH/Mz) |22j . 

To emphasize the main features of the quantum corrections to the neutrino mass 
matrix, we will concentrate in what follows on a see-saw model with just one right- 
handed neutrino with mass Mmaj- Then, with this assumption, the neutrino Yukawa 
couplings Fjf are 3- vectors. In this scenario the neutrino mass matrix at the scale of the 
right-handed neutrino mass Mmaj is given by: 

ylylT 

[^''r'' = ir^ , (16) 

which has rank 1 and thus only one non- vanishing eigenvalue. 

On the other hand, quantum effects introduce corrections to the neutrino mass matrix 
yielding at low energies K^^{mH) = [n^^Y^^'^ + 6k,^^, where the correction reads, in the 
leading-log approximation, 

S^'' = --^/3.n(M,,aj)l0g^^ . (17) 

Ibn^ niH 

Using the explicit form of the /3 function in the Appendix, it follows that this correction 
can be schematically written as: 

6k^' ~ BiaK''^ + K^'^Bj^ + bK^^ . (18) 

Here Bia denote 3x3 matrices whereas 6 is a number. The first two terms generalize the 
well known correction to the neutrino mass matrix in the Standard Model including the 
dimension-5 Weinberg operator. However, the term proportional to k^^ does not have 
any correspondence in the Standard Model and, as we will see, introduces new qualitative 
features. The coefficient b explicitly reads: 

' ,2A5log^^, (19) 



IGtt^ rriH 



which depends hnearly on the coefficient of the potential term A5(<l>j$2)(^i^2), while 
only logarithmically on the ratio between the scale of the right-handed neutrino and the 
overall scale of the extra scalars H^, H^, A^. 

The neutrino mass matrix Ai^ = [AiuY'^'^^ + ^-Mi, can be diagonalized using perturba- 
tion theory, giving as a result the eigenvalues rrii = m\ + Snii with m\ the eigenvalue 
at tree level and 6mi the ffist order correction. 

At lowest order in perturbation theory, taking into account only the tree level mass 
term, there is only one non-vanishing neutrino mass eigenvalue: 



mr = W^K'\''- (20) 



,(0) _ ^ 1^^112 



/iWmaj 



On the other hand, the third column of the leptonic mixing matrix reads; 

4°^ = ^ ' (21) 

\ u \ 

while the first two columns are undefined, due to the degeneracy of the corresponding 
neutrino mass eigenvalues. In this expression, \Yj^\ = (Yli I^iaI'^Y^'^- 

The correction to the neutrino mass eigenvalues due to the perturbation 6k^^ is given 
by: 

<5m, = ^Re[(f/(°)^5«:i^f/W)..], (22) 

which slightly modifies the value of the heaviest neutrino mass eigenvalue: 

^2n2" 



(5mq = — -— — Re 



7;2 r ryitya 



2M^aj 



Y,} 



(23) 



More importantly, this correction is also non- vanishing for 6m2, thus generating ra- 
diatively a second neutrino mass eigenvalue. This is in contrast to the widely studied 
case of the Standard Model extended with a single right-handed neutrino, where there 
is only one non-vanishing neutrino mass eigenvalue, even after taking into account the 
renormalization group running. \j 

To show this, we write explicitly the radiative correction to the next-to-lightest neu- 
trino mass eigenvalue: 



Sm, = ^-^Re[UJS6^lMt'], (24) 



which crucially depends on the second column of the zero-th order leptonic mixing ma- 
trix. Since the matrix f/^^^ is unitary, the vector f/^2 should satisfy ^qUJ^2Uq3 — 0' 
So t/g2^g2 = 1- A vector that satisfies those properties can be easily constructed from 



""^ There are, however, tiny finite corrections arising from two-loop diagrams involving W bosons 



the vectors Y^ and Y^ using the Gram-Schmidt process. Starting with U-^ given by 
eq. (I2T]) one finds: 



U, 
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I arg(-A5) 



(25) 



where 



No 



Y^'fy^ 
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1/2 



(26) 



Substituting into the expression for 61712 we find 
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M, 
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m^ 



(27) 



(Note that the phase in eq. fl25|) has been chosen to yield 7712 real and positive.) It is 
apparent from this expression that in order to generate a non-vanishing neutrino mass 
eigenvalue it is necessary the misalignment between the Yukawa couplings Y,} and Y^, or 
in more physical terms, it is necessary the existence of new sources of fiavour violation 
in the neutrino sector. These new sources necessarily generate, through quantum cor- 
rections, off-diagonal elements in the charged lepton Yukawa coupling Y^, which in turn 
induce a contribution to the lepton flavour violating processes. Nevertheless, as explained 
in Section [2] this contribution is suppressed by the large mass of the extra Higgs particles, 
and can be consistent with experiments if the extra particles are sufficiently heavy. 

Furthermore, it is interesting to note that, under some well motivated assumptions, 
the hierarchy between the tree level mass m^ and the radiatively generated neutrino mass 



m2 can be fairly mild. For instance, taking the typical values IA5I ~ 1, Mj 



maj 



and rriH ~ 1 Te V and assuming non-aligned neutrino Yukawa couplings with | Y^^ 
one obtains for the ratio between the two heaviest neutrino mass eigenvalues 
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(28) 



which yields a mild mass hierarchy, in qualitative agreement with the experimental data. 
Note that, whereas the overall scale of the light neutrino masses depends linearly on the 
inverse of the heavy right-handed neutrino mass, the ratio between the two heaviest neu- 
trino mass eigenvalues depends only logarithmically with the masses of the new particles. 
As a consequence, the result in eq. ( l28l) is fairly insensitive to the exact values of the 
masses of the heavy particles. 

In the previous analysis we have assumed for simplicity that only one right-handed 
neutrino participates in the neutrino mass generation. In the more realistic case where 
there are several right-handed neutrinos, the tree level contributions to all neutrino mass 
eigenvalues will be non- vanishing. Nevertheless, as discussed in the Introduction, if the 



9 



neutrino Yukawa coupling Y^ has hierarchical eigenvalues, as suggested by the observed 
hierarchies in the quark and charged lepton masses, then the neutrino mass hierarchy 
generated (at tree level) by the see-saw mechanism is in general several orders of magni- 
tude larger than the one inferred from experiments. Therefore, the radiatively generated 
contribution to the next-to-lightest neutrino masses by the presence of the second Higgs 
doublet will dominate over the tree level contribution, and the conclusions presented 
above will still hold. 

In an extended scenario with Nh Higgs doublets and one right-handed neutrino, the 
radiatively induced next-to-lightest neutrino mass receives A^^x {Nh — 1)/2 contributions, 
each of them proportional to the coefficient of the term A5*($jj$i)($|^$i), a,b = 2...Nh, 
in the Higgs potential. Therefore, in this case the radiatively generated neutrino mass is 
enhanced. 

It is amusing to speculate that adding more Higgs doublets to the particle content of 
the model may also be relevant to understand the observed pattern of neutrino mixing 
angles. In a model with Nh Higgs doublets, only the Higgs that acquires a vacuum 
expectation value, $1, will contribute to the tree level mass. Assuming that this is 
the largest mass, ms = \/Aml^^, it follows from eq. ( 12T|) that Ui^ oc Y^*. Therefore, 
if there is any pattern in the neutrino Yukawa coupling yjj, stemming e.g. from an 
underlying flavour symmetry, then this pattern will be inherited by t/js, thus providing 
an explanation to the apparent structure of the last column of the leptonic mixing angle: 
|f/i3| ~ 0, If/23 1 — If/as]. In contrast, there are N^ — 1 Higgses which contribute via 
quantum effects to the generation of the solar neutrino mass scale, m2 = y^Am^ and of 
the second column of the leptonic mixing matrix, Ui2- As a consequence, even if there is 
a structure in each of the neutrino Yukawa couplings, YJ^, a = 2, ...Nh, the generated f/j2 
will be structureless, since it receives contributions from all these Yukawa couplings. This 
is in rough agreement with observations, which reveal that the three entries in Ui2 are 
all (9(0.1), without displaying any remarkable structure (or, alternatively, the solar angle 
is neither maximal nor zero). Therefore, in the Standard Model extended with right- 
handed neutrinos and several Higgs doublets, the last column of the leptonic mixing 
matrix is expected to display a "hierarchical" structure, whereas the second column, an 
"anarchical" structure [23], in qualitative agreement with the data. 



4 Comparison to the two right-handed neutrino 
model 

The scenario discussed in this paper leads to a dimension-5 operator which is identical 
to the one generated by the Standard Model (with a single Higgs doublet) extended by 
two heavy right-handed neutrinos. There are however some conceptual differences in the 
way these two scenarios reproduce the observed neutrino data, which we discuss here. 

Let us first demonstrate the equivalence of the two Higgs doublet model extended 
with one right-handed neutrino (2HD-1RHN model) and the Standard model extended 
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with two right-handed neutrinos (1HD-2RHN model). At low energies, the coefficient of 
the dimension five operator generated in the 2HD-1RHN model reads, following eqs. flT3|) . 

K"(m^) ^ -^ fcy;^ + B,,Y,^Y,'^ + YX'^Bl + bY^vf] , (29) 

-''^maj L 

which can be recast as 

t^\mH) ^ j^ \{Y,' + B,,Y:){Y^' + B,,Y,^f - B.^Y^Yf B^^, + bY^Yf] . (30) 

Neglecting the term of 0{B^) and defining Y^} = Y^ + BiaY^, one obtains the following 
low energy neutrino mass matrix 

, yifi^ bY,^Yf\ v^ , , 

-'•'-'inaj -'"maj 

This result is formally identical to the effective operator which arises in the low en- 
ergy limit of the 1HD-2RHN model. More specifically, in this model the high energy 
Lagrangian reads, in the basis where the 2x2 right-handed neutrino mass matrix is 
diagonal, 

- Z:'^-'^™ = {Yu)4Li^R,^ - \mi1^ri^ri - ^M^uguR^ + h.c. , (32) 

which leads, when Mi, M2 ^ f , to the effective neutrino mass matrix 

■^^ - V Ml + M2 J 2 ' ^^^^ 

being Yi^2 column vectors defined as Yi = (Yu)ii, Y2 = {Yi^)i2. By comparing eqs. (l3T]l 
and ( 133|) it follows that, from the point of view of the neutrino mass generation, the 2HD- 
IRHN model is equivalent to the 1HD-2RHN model with the following correspondence 
among parameters 

{Fi, Y2, Ml, M2} ^ {Y,\ Y^, M^aj, M^ajM . (34) 

This correspondence allows to write explicit expressions for the most general Yukawa 
couplings Y^, Y^ which lead to the neutrino masses m2^m^ and the leptonic mixing 
matrix U. Using the results of [25], one easily finds: 

Yl = ^^/M~{V^2Coseu:, ± ^/W^smeU*,) , (35) 

^, = ^y^(_y^sin^V;^ ± ^W.coseu:,) , (36) 

11 



where 6* is a complex angle which parametrizes the family of Yukawa couplings compatible 
with the low energy neutrino data (note that this parametrization may fail if the term 
of 0{B'^) in eq. (150]) can not be neglected). Finally, the Yukawa coupling with the Higgs 
$1 is Y^ = (1 - fin)n' - B^^Y^. 

Furthermore, this correspondence allows to better appreciate the advantages of the 
2HD-1RHN model over the 1HD-2RHN model in the generation of a mild mass hierarchy. 
In the 1HD-2RHN model, the neutrino mass hierarchy is essentially given by 



mg |F2rMi 
m2 |Fi|2M2 



(37) 



therefore the only possibilities to generate a mild neutrino mass hierarchy are i) 1^2^ ~ 
|Fi|2, M2 ~ Ml, or ii) IF2I > l^il with M2/M1 ~ |Yi|Vl^2p- In view of the observed 
large hierarchies in the quark and charged lepton Yukawa eigenvalues, in a model with 
two right-handed neutrinos one expects \Y2\ 2> \Yi\, which hence requires a huge hi- 
erarchy between the two right-handed neutrino masses in order to render a mild light 
neutrino mass hierarchy. More concretely, if the neutrino Yukawa eigenvalues have a 
similar hierarchy as the up-type quark masses, |l2|/|^i| ~ mt/nic ~ 150, then it is re- 
quired M2/M1 ~ 20000. On the other hand, if the hierarchy is similar to the down-type 
quark masses, |l2|/|^i| ~ fnb/rris ~ 40, then it is required M2/M1 ~ 1600. Whereas such 
large hierarchies in the right-handed neutrino masses cannot be precluded, it is difficult 
to conceive that in the decoupling limit the large hierarchies in the Yukawa couplings 
cancel almost exactly a huge hierarchy in the right-handed neutrino masses to generate 
at low energies a light neutrino mass hierarchy of ~ 6, as inferred from experiments. 

This drawback is very naturally circumvented by the 2HD-1RHN model which as we 
argued above is equivalent, concerning the neutrino mass generation, to the 1HD-2RHN 
model. In the equivalent 1HD-2RHN model, the Yukawa couplings Yi and Y2 are naturally 
of the same size, since in the original 2HD-1RHN they correspond to Yukawa couplings 
to the same generation of right-handed neutrinos. Furthermore, in the corresponding two 
right-handed neutrino model, the masses Mi and M2 naturally present a mild hierarchy, 
given by the factor 1/b ~ 0{1 — 10). As a result, the 2HD-1RHN is equivalent to a 
1HD-2RIIN model which naturally fulfills the conditions i) to reproduce the observed 
mild neutrino mass hierarchy. 

Another important difference between the 1HD-2RHN model and the 2HD-1RHN 
model concerns the possibility of observing other phenomena at low energies apart from 
neutrino masses. It is well known that, in its simplest version, the 1HD-2RHN model does 
not have any other observable low energy prediction apart from the tininess of neutrino 
masses. In this model the scale of lepton flavour and lepton number violation both coin- 
cide with the scale of the right-handed Majorana neutrino masses, which are postulated 
to be much larger than the electroweak symmetry breaking scale. As a consequence, 
the rates of all flavour and lepton number violating processes are inversely proportional 
to the heavy right-handed neutrino mass resulting in tiny rateso In contrast, in the 



^A notable exception is the 1HD-2RHN scenario where the right-handed neutrinos form a pseudo- 
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2HD-1RHN model, apart from the lepton flavour associated to the right-handed neu- 
trino couphngs there exist another source of lepton flavour violation associated to the 
charged lepton couplings to the second Higgs doublet, inducing rates for the rare lepton 
decays suppressed by the heavy Higgs masses. If the additional scalar degrees of freedom 
have masses not far from the electroweak symmetry breaking scale, the induced rates of 
the lepton flavour violating processes could be large enough to be observed in experi- 
ments. A more detailed discussion about the prospects to observe the process /i — )■ 67 in 
experiments will be presented in Section El 

5 Corrections to the mixing angles and discussion of 
sin ^13 

Below the right-handed neutrino mass scale, the neutrino mixing angles receive radiative 
corrections with two different origins. First, the change in the entries of the neutrino 
mass matrix due to the RGE running, eq. ( !T8|) . generates a correction to the leptonic 
mixing matrix given by: 



where 



5U^ = U^^^T , (3^ 






T^ = ^p ^p if ^7^ J- (40) 



However, this is not the physical leptonic mixing matrix measured by experiments, 
since the RGE running also modifles the structure of the charged lepton Yukawa cou- 
plings. More speciflcally, if the charged lepton Yukawa coupling Y^ is diagonal at Mmaj, 
the radiative corrections induced by Y^ will generate at low energies off-diagonal entries 
in Y^. It is then necessary to redeflne the charged lepton flelds in order to render a 
diagonal charged lepton Yukawa coupling, namely II -^ Ve^L, ^r — ^ V^cr, where V^, 
V^ follow from the singular value decomposition, Y} = V/diag(?/gi,?/g2,Z/e3)K^^- This 
redeflnition introduces an additional correction to the leptonic mixing matrix given by 

SUy^ = (K^ - l)^t/(°) . (41) 



Dirac pair with masses 0(100 — 1000) GeV. In this case, the Yukawa couphngs can be sizable while 
correctly reproducing the tininess of the neutrino masses. As a consequence, the rates for ^ — )- 67 and 
neutrinoless double beta decay can be largely enhanced, possibly allowing their observation in the next 
round of experiments [55] . 
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The matrix V^^ can be explicitly calculated from the /3-functions of the charged lepton 
Yukawa couplings. Using 



ivitl 



Y'Y 



\mH 



Y}Y}^ - i(%i;^^ + Y^^U log ^^ 



167r2 



we obtain 



(K 



1 (%n^^ + nV^k- M_ 



e )ij 



167r2 (yi)^ - (2/e\) 



M2 



log- 






m^ 



^7^ J • 



(42) 



(43) 



Therefore, summing up the two contributions, the leptonic mixing matrix at low 
energies reads, in the physical basis. 



[/(I) = vyU^""^ + [/(°)T 



(44) 



We are particularly interested in the correction to the last column of the leptonic 
mixing matrix, which in general yields a non- vanishing contribution to sin ^13 and a 
deviation from the maximal atmospheric mixing which may be observed in experiments. 
Concretely, the correction to the third column of the leptonic mixing matrix due to the 
running of k reads 



(5f/.).3 = (f/^°^T) 



Y, 



1* r 



i3 



ReiY^^B^aY^) , Um(6*(F,2ty^i) 
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(45) 



while the contribution from the rediagonalization of the charged lepton Yukawa coupling 
reads: 
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Yji3 



IGVT' 
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(46) 



A quantity of particular interest is the angle ^13, which is constrained by present 
experiments to be small. It is interesting that radiative corrections can generate in this 
model a fairly large value of ^13, possibly at the reach of the planned experiments, even 
if its tree-level value vanishes. Summing up the contributions from eqs. (HSj) and 046 p . 
and neglecting terms cubic in the charged lepton Yukawa couplings, we obtain that the 
radiatively induced value of [/13 is 



2* 



5U 



13 
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(47) 
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which is, as the ratio m2/m^, suppressed by the loop factor but enhanced by the large 
logarithm of the ratio of the Majorana mass over the heavy Higgs mass. As a result, the 
radiatively generated 6is can be as large as ~ 0.2 if any of the entries in the bracket is 
~0(1). 

So far we have considered only the corrections to the neutrino mixing angles from 
the running below the right-handed neutrino mass scale. However, in a large class of 
models the cut-off of the theory lies at higher energies and additional contributions to 
the leptonic mixing may arise from radiative corrections between the cut-off scale A and 
the right-handed neutrino mass scale M^aj- The corrections for this case can be derived 
using eqs. fH5l) . P6|) and the substitution rules given in Appendix I A. 2 1 the result being: 
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13 
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(48) 



Quantum effects also induce corrections to the atmospheric mixing angle, leading to 
deviations to the maximal mixing even if ^23 = "^ /^ at tree level. It is interesting that 
if the neutrino Yukawa couplings are the dominant source of flavour violation in the 
leptonic sector, then a correlation arises between the deviations of t/23/^33 and ?7i3 from 
their corresponding values at the cut-off scale. 

In this limit, the radiative corrections to the last column of the leptonic mixing matrix 
are dominated by the first line of eq. P7|) . which can be schematically written as: 



U. 



i3 



r{0) 



ez)U^' + e2f/,2 
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(49) 



It can be checked that to first order the ratio U23/U33 does not depend on €3. Then, using 
the equation for U13 to eliminate €2 it follows that: 
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(50) 



Concretely, in the case when at the cut-off scale the atmospheric mixing angle is 
exactly maximal and 613 vanishes, at low energies the elements of the leptonic mixing 
matrix approximately satisfy 

^ - 1 ~ 2v^[/i3 , (51) 

(^33 

which can be recast as 



tan6'23 ~ |l + 2V2sin6' 



136 



or 



^23- - ^ 72 sin ^13 cos (5. (52) 



15 



CO 

o 



0.15 




0.05 



-0.15 -0.1 -0.05 0.05 0.1 0.15 



'<2 

o 



0.15 


1 1 


1 


1 1 


0.1 


^^4.. 




, ^el^^^ 




"•^^^^fem^ 4-^ . 


',' -, .." 


^M^^ 


0.05 


~'.^^^^nH^ 


i|4l|ij 


y^"^ 




^'''^My^m 


^>" 


y^ ^ 












' ' '-. ' '?*^^ J 


l^s. -/ 




n 


' ,!•■.«&«- 


^^^n^^BCTr r- 


i > 



-0.15 -0.1 -0.05 0.05 0.1 0.15 

^23 - T 



Figure 1: Scatter plots showing | sin 6*13 cos 5| against 6*23 — ^ at low energies for ran- 
dom choices of high energy parameters consistent with the measured neutrino oscillation 
parameters. We have assumed tri-bi-maximal mixing at the cut-off scale, being the devi- 
ation from 6*23 = vr/4 and ^13 = at low energies only due to the radiative corrections as 
described in the main text. The cut-off scale is A = Mmaj = 10^^ GeV in the left panel 
and A = 10^^ GeV in the right panel. 



If there are additional sources of lepton flavour violation, then the low energy values 
of 623 and ^13 are expected to deviate from this relation. This is illustrated in the scatter 
plots shown in fig. [H which have been obtained by the numerical one loop integration of 
the RGBs of the 2HDM extended by one right-handed neutrino. We assume in the plot 
rriH = 3 TeV, Mmaj = lO^'^ GeV and tri-bi-maximal mixing at a cut-off scale, which we 
take A = M^aj (A = 10^^ GeV) in the left (right) panel. Quantum effects generate a 
non- vanishing value for ^13 and ^23 — f, mostly due to the RGE effects of A5 and Ag, as 
follows from eqs. p7|) and pHjl . In the plot we fix IA5I = 0.5 and we take random values 
with jAgl < 0.45, in order to preserve the perturbativity of the quartic couplings in the 
renormalization group running. To investigate the impact of the charged lepton mixing 
in the correlation we have adopted the ansatz Y^ = VY^, where l^ is a general unitary 
matrix with random angles and phases. Furthermore, since the effect of the charged 
lepton Yukawa couplings on the corrections to the leptonic mixing matrix is proportional 
to Tr(FjF„^"l'), we have taken in the scatter plot 1(^^^)33/(^^)33! < 0.05 (red points), 
l(>;')33/(>;')33| < 0.15 (green points) and \{Y^)33/iY^)33\ < 0.3 (blue points); the effects 
of the down quark Yukawa couplings have been neglected in this analysis, although their 
role is completely analogous. All the points in the plot reproduce the neutrino oscillation 
parameters within their experimental errors. 

It is apparent from the plots that when the charged lepton Yukawa couplings have a 
negligible effect on the running (corresponding to |(^„^)33/(yu^)33| ^ 1), there is a fairly 
strong correlation between the radiatively generated 6*13 and 6'23 — ^. When the cut-off is 
A = Mmaj, the numerical results are in good agreement with eq. f lS2]) . shown as a black 
solid line in the plot. In contrast, when A = 10^* GeV there is a larger spread of the 
points, due to the additional RGE effects between A and M^aj- Besides, in this case the 
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numerical results do not agree with eq. fl52|) . since the RGE running between A and Mmaj 
generates a non- vanishing (and negative) shift of U23/U33 at the scale M^aj, even if U13 
still vanishes. This produces, following eq. fISUl) . the shift of the points to the left of the 
black solid line. 

To summarize, from our analytical and numerical analysis, it follows that in the 
2HDM extended with one right-handed neutrino it is generally expected a deviation of 
the atmospheric angle from the maximal value which is comparable to the reactor angle 

TT 






^13 . (53) 

unless the CP violating phase S is very close to 7r/2. 



6 Lepton flavour violation 

In the general 2HDM extended with RH neutrinos, one generically expects a misalign- 
ment in the charged lepton Yukawa couplings which will lead to new phenomena at low 
energies, apart from neutrino masses, in contrast to the standard see-saw scenario with 
just one Higgs doublet. This misalignment will generically arise already at tree level. 
However, even if the charged lepton Yukawa couplings are aligned at the cut-off scale A, 
radiative corrections from the neutrino Yukawa couplings from the RGE running above 
the Majorana mass scale will introduce off-diagonal entries in both charged lepton Yukawa 
matrices. Note that the radiative generation of neutrino masses requires a misalignment 
in the neutrino Yukawa couplings, hence some amount of flavour violation is necessarily 
generated via quantum corrections in the charged lepton sector. 

To calculate the minimum amount of lepton flavour violation in the charged lepton 
sector, we assume that Y^, Y^ are diagonal at the cut-off scale A > Mmaj- Then, due 
to the radiative corrections from the neutrino Yukawa couplings Yj^, Y^, both charged 
lepton Yukawa couplings become non-diagonal at the scale Mmaj- As discussed in the 
previous section, we now redefine the charged lepton fields in order to bring the Yukawa 
coupling Y^ into its diagonal form. As a result, the off-diagonal elements of the charged 
lepton Yukawa coupling read, at the Majorana mass scale: 

log- 



Below the Majorana mass scale the charged lepton Yukawa couplings are also affected by 
the quantum effects, however the off-diagonal elements at low energies are still given by 
the previous expression, up to second order effects. 

To estimate this contribution we assume Y^ = C,eY^, with at the cut-off scale. With 
this ansatz, eq. ( I54l) reads: 



log- 



A 



(K^^i;')l2L/ . = -^T^iYu + ^eY^-^eY^'* + If *)2i;22 • (55) 

17 



Inserting this contribution into eq. (|8]) one obtains the approximate lower bound: 



2 



BR(/. -> e7) > ^ [^^ 1 K\ + ^eY^W'Kl - CX.\' 



mf\ lAfilw 



mlJ mjj 



(56) 



which is saturated when the charged lepton Yukawa couphngs are ahgned and when the 
Yukawa couphngs of the heavy Higgs to the quarks are neghgible. 

To estimate the size of this lower bound, we will assume neutrino Yukawa couplings 
maximally misaligned with the form Y^} = -^(0,1,1)"^, Y^ = -^(1,1,-1)^, being yi 
and 2/2 the corresponding norms. With this choice, we obtain, in the limit |,^e| ^ 1 and 
niH > Mz and taking A = lO^Mmaj, 

BR(/i -> 67) > 3 X 10-1^ X |A62/i?/2CT {j^y ■ (57) 

Given that this bound is very conservative, the observation of the process fi -^ e'-f may be 
at the reach of the MEG experiment, which aims to BR(/i — )■ 67) > 10~^^ |27], provided 
the couplings are sizable and provided the extra scalar degrees of freedom are not too 
heavy. 

7 Conclusions 

We have considered in this paper an extension of the Standard Model by one extra Higgs 
doublet and one or more Majorana right-handed neutrinos, including in the Lagrangian all 
terms compatible with the Standard Model gauge symmetry. We have calculated, using 
a renormalization group approach, the quantum corrections to the neutrino parameters 
under the assumption that the right-handed Majorana mass scale is much larger than the 
mass of the Higgs mass eigenstates. We have argued that if the neutrino Yukawa couplings 
are misaligned, the radiatively generated contribution to the mass of the next-to-heaviest 
neutrino can be much larger than the tree level mass. Furthermore, for reasonable choices 
of the parameters of the model, the radiatively generated mass of the next-to-heaviest 
neutrino is a factor of a few smaller than the mass of the heaviest neutrino. Since the mass 
hierarchy depends only logarithmically on the masses of the extra degrees of freedom, 
this conclusion is fairly insensitive to the scales at which the new physics appears. 

Hence, in this simple model two puzzles in neutrino physics can be simultaneously 
explained. First, the smallness of the neutrino masses is explained by the see-saw mecha- 
nism. Secondly, the mild hierarchy between the atmospheric and the solar neutrino mass 
scales is explained by the radiative origin of the mass of the next-to-heaviest neutrino, 
which is suppressed by the loop factor but enhanced by the large logarithm of the ratio 
between the heavy Majorana mass scale and the heavy Higgs scale. Furthermore, by 
making the heavy Higgs scale sufficiently large, all the successes of the Standard Model 
can be preserved, since all low energy effects of the extended Higgs sector are suppressed 
at least by two powers of the heavy Higgs mass. 



The misalignment in the Yukawa couphngs, necessary for the radiative generation of 
the solar neutrino mass scale, amounts to new sources of lepton flavour violation which 
also modify the structure of the leptonic mixing matrix through the renormalization group 
evolution. Therefore, we expect in this model deviations from the maximal atmospheric 
mixing and from a vanishing ^13 due to quantum effects. We have carefully calculated 
these corrections and we have found that the radiatively generated angle 9i^ can be large 
enough to be measured in present and future experiments. 
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A Appendix 



A.l Quantum corrections below Mj 



maj 



The one-loop /3 functions of the multi-Higgs doublet model, including the dimension five 
operator which yields neutrino masses, have been derived in [2B]. The fi functions of 
the charged lepton Yukawa couplings, Y^, and the dimension five operators k°-^ read, for 
energy scales below the right-handed neutrino Majorana mass scale, 



/5ya 



9 2 



^^" ) y: + 



3Tr ( F:%^ + yj^F, 



ct 



Tr {y:y:^) 



Y' 



+ y:y:^y: + ^:Y:^Y:, 



{51 



(3.. 



ab — 



y^'y^^^k"'' + K^^Y^^Y^^^)^ +2 y;i;''tfi:"" + K^* (1;%'^^) 



vT 



Y^'Yf{n'' + 



K 



bc\ 



[K,- 



+ o(>;%'0' 



+ 



+ K 



3Tt{Y:y:^ + YfY,^) + Tr(r;ty; 
STriY^Y:^ + Y^^Y^) + Tr(F,^ty^ 



.cfe 



- Sg'K''' + 2A,,M/€' 



cd 



(59) 



where summation over repeated indices is understood and the quartic couplings A are 
defined by \/ D |A,fe,,(<Dt$6)($t$rf). 

In the case of only one right-handed neutrino and two Higgs doublets, n^^^rriH) can 
be approximately written at the leading log at any mass scale rriH < Mmaj in the form 
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of eq. f lTSj) which we repeat here for completeness 
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where we have defined flavour matrices i?ii, B12 and the complex number h by 



167r2 






■i?ii = -in'n'^ + lYeYe' - 3Tr(rjFjt + y^V^i) _ Tr(i;ity^i 



167r2 



logHsiM 



B,, = 2YX'^ - 3Tr(yJyJt + Y^%') - Tt{Y,'X') - 2\ 



167r2 






-2A, . 



2^ 



Ai 



(62) 
(63) 

(64) 



A. 2 Quantum corrections above M, 



maj 



In the case that the cut-off scale of the theory, A, is larger than Mmaj, the relevant 
matricial couplings of the leptonic Lagrangian are the charged lepton Yukawa couplings, 
Y^, the neutrino Yukawa couplings, Y^, and the right-handed Majorana mass matrix 
Mm- The corresponding /3-functions are: 
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The running above Mmaj modifies some of the expressions that we have derived in this 
paper. The effects of the running can be easily incorporated in our results by substituting 



a\rbT 



Bla — > Bia + i?i„ , 

More concretely, the values for B^^ read: 
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Note that the running above Mmaj does not modify the value of b nor m2, cf. eq. (|2 

Besides, the running above Mmaj also affects the structure of V^^ at low energies. This 
can be taken into account with the following substitution in the relevant formulas: 

(3y. log ^^ ^ I3y. log ^^ + /3^. log -^ . (71) 

References 

[1] For a review, see M. C. Gonzalez- Garcia, M. Maltoni, Phys. Rept. 460 (2008) 1-129. 

[2] M. Mahoni, T. Schwetz, PoS IDM2008 (2008) 072. 

[3] P. Minkowski, Phys. Lett. B 67 (1977) 421; M. Gell-Mann, P. Ramond and R. 
Slansky, Proceedings of the Supergravity Stony Brook Workshop, New York 1979, eds. 
P. Van Nieuwenhuizen and D. Freedman; T. Yanagida, Proceedinds of the Workshop 
on Unified Theories and Baryon Number in the Universe, Tsukuba, Japan 1979, eds. 
A. Sawada and A. Sugamoto; R. N. Mohapatra, G. Senjanovic, Phys. Rev. Lett. 44 
(1980)912, ibid. Phys. Rev. D23 (1981) 165; S. L. Glashow, The Future Of Elementary 
Particle Physics, In *Cargese 1979, Proceedings, Quarks and Leptons*, 687-713 and 
Harvard Univ. Cambridge - HUTP-79-A059 (79,REC.DEC.) 40p. J. Schechter and 
J. W. F. Valle, Phys. Rev. D 22 (1980) 2227. 

[4] J. A. Casas, A. Ibarra and F. Jimenez- Alburquerque, JHEP 0704 (2007) 064. 

[5] E. Ma, Phys. Rev. D73, 077301 (2006). 

[6] See, for example, J. Kubo, E. Ma, D. Suematsu, Phys. Lett. B642 (2006) 18-23; 
C. Boehm, Y. Farzan, T. Hambye, S. Palomares-Ruiz, S. Pascoli, Phys. Rev. D77 
(2008) 043516; M. Aoki, S. Kanemura, O. Seto, Phys. Rev. Lett. 102 (2009) 051805, 
Phys. Rev. D80 (2009) 033007; J. Kubo, D. Suematsu, Phys. Lett. B643 (2006) 
336-341; N. Haba, K. Tsumura, JHEP 1106 (2011) 068. 

[7] W. Grimus, H. Neufeld, Nucl. Phys. B325 (1989) 18. W. Grimus and L. Lavoura, 
Phys. Lett. B 546 (2002) 86. 

[8] W. Grimus, H. Neufeld, Phys. Lett. B486 (2000) 385-390. 

[9] T. D. Lee, Phys. Rev. D8, 1226-1239 (1973). 

[10] For a recent review, see G. C. Branco, P. M. Ferreira, L. Lavoura, M. N. Rebelo, 
M. Sher, J. P. Silva, |arXiv: 1106.0034, [hep-ph]]. 

[11] M. E. Peskin, T. Takeuchi, Phys. Rev. Lett. 65 (1990) 964-967. 

[12] J. F. Gunion, H. E. Haber, Phys. Rev. D67 (2003) 075019. 

21 



[13] H. E. Haber, D. O'Neil, Phys. Rev. D83 (2011) 055017. 

[14] G. D'Ambrosio, G. F. Giudice, G. Isidori and A. Strumia, Nucl. Phys. B 645, 
155 (2002); A. Pich and P. Tuzon, Phys. Rev. D 80, 091702 (2009); A. J. Buras, 
M. V. Carlucci, S. Gori and G. Isidori, JHEP 1010, 009 (2010); C. B. Braeuninger, 
A. Ibarra and C. Simonetto, Phys. Lett. B 692 (2010) 189; M. Jung, A. Pich and 
P. Tuzon, JHEP 1011 (2010) 003, Phys. Rev. D 83 (2011) 074011. 

[15] M. L. Brooks et al. [MEGA Collaboration], Phys. Rev. Lett. 83 (1999) 1521. 

[16] P. Paradisi, JHEP 0602 (2006) 050. 

[17] J. Hisano, S. Sugiyama, M. Yamanaka et al, Phys. Lett. B694 (2011) 380-385. 

[18] J. Laiho, E. Lunghi and R. S. Van de Water, Phys. Rev. D 81 (2010) 034503. 

[19] K. Nakamura et al. [ Particle Data Group Collaboration ], J. Phys. G G37 (2010) 
075021. 

[20] A. J. Buras, S. Jager, J. Urban, Nucl. Phys. B605 (2001) 600-624. 

[21] E. Lunghi and A. Soni, JHEP 0709 (2007) 053. See also A. Lenz, U. Nierste, JHEP 
0706 (2007) 072. 

[22] J. A. Casas, J. R. Espinosa, A. Ibarra, I. Navarro, Nucl. Phys. B573 (2000) 652- 
684, Nucl. Phys. B569 (2000) 82-106; J. R. EUis, S. Lola, Phys. Lett. B458 (1999) 
310-321; P. H. Chankowski, W. Krolikowski, S. Pokorski, Phys. Lett. B473 (2000) 
109-117; S. Antusch, J. Kersten, M. Lindner, M. Ratz, Nucl. Phys. B674 (2003) 
401-433. 

[23] S. T. Petcov, S. T. Toshev, Phys. Lett. B143 (1984) 175; K. S. Babu, E. Ma, Phys. 
Rev. Lett. 61 (1988) 674; E. Ma, Phys. Lett. B456 (1999) 48-53. 

[24] L. J. Hall, H. Murayama, N. Weiner, Phys. Rev. Lett. 84 (2000) 2572-2575. 

[25] J. A. Casas and A. Ibarra, Nucl. Phys. B 618 (2001) 171. A. Ibarra, G. G. Ross, 
Phys. Lett. B591 (2004) 285-296. 

[26] A. Ibarra, E. Molinaro, S. T. Petcov, |arXiv:1103.62T7l [hep-ph]]. 

[27] A. Maki, AIP Conf. Proc. 981 (2008) 363-365. 

[28] W. Grimus, L. Lavoura, Eur. Phys. J. C39 (2005) 219-227. 



22 



